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ABSTRACT 


Historically, peristalsis, the rhythmic, traveling wave like 
pumping of bodily organs such as the ureter, has been mechanically 
Simulated by the extended longitudinal motion of rollers squeezing 
fluid along a flow channel. Blood pumps are an example. The 
experimental work associated with this thesis has developed a 
peristaltic pump wherein pumping is produced by the limited 
transverse motion of electromechanical transducers, The transducers 
are arranged in a line and excitation is time phased so that 
their sequential transverse motions produce a traveling deflection 
wave along the pump flow tube, 

The pump is a simple device made from transformer cores, 

Bach transformer serves as a transducer; the "E", with windings 
on the center leg, provides AC and DC excitation; the "I" serves 
as the actuator which squeezes the pump tube and produces flow, 
The magnetic action is much like a solenoid in that the "E" and 
"I" serve as a magnetic circuit and the magnetic field in the air 
gaps produces the tube squeezing force, 

The theory developed in the thesis involves the electromechanical 
characteristics of each pump transducer and the fluid mechanical 
flow characteristics for the pump. The final tying together of 
theory and experiment is a comparison of theoretical pump pressure- 
flow characteristic curves and experimental pump operating points, 
Operating points do conform fairly well to the pressure-flow 
characteristics, but insufficient operating points were determined 
to state that the pressure-flow curves accurately predict pump 
operation, 

Production of peristaltic pumping by transverse deflections 
has great potential value. If a peristaltic pump is ever to 
be used as a replacement for an organ of the body, longitudinal 
rollers are not practical. The transverse deflection pump is 
simple and also might be miniaturized. In place of transformers, 
piezoelectric crystals could be vsed as transducers. With a 
DC voltage excitation, crystals could produce great forces over 
short displacements, Miniaturization and improvements as recommended 
in the thesis are worthwhile areas for future experimentation 
and development, 
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Name and Unit of the Symbol 
Transducer outer leg einVelnees: meters 
Transducer exciter center leg thickness, meters 
fometic flux density, webers/m* 
Phase velocity, m/sec 
Transducer exci.ter width, meters 
Damping coerficient 
Vector differential operator 
Frequency, hertz 
Electric force, newtons 
Fluid pressure force, newtons 
Spring force, newtons 
Magnetic flux, webers 
Gravity acceleration, fi aeee 
Normalized channel. width, r/R 
Magnetic field intensity, amp-turns/m 
Pump suction head, meters 
Pump discharge head, meters 
Time varying current, amps 
Steady state current, amps 

ai 
i)% 
Current Density, Pee 
Wave number, 211/A 
Springs constant, newtons /1m 
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AP 
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Name and Unit of the Symbol 





Petuacor width, d/ 11, meters 
Inductance, henries 

Wavelength, meters 

Actuator mass, kg 

Dynamic viscosity, merdonseee | 
Permitivity of free space, henries/m 
Exciter winding turns, turns 
Kinematic viscosity, A] p 

Fluid pressure, neumeen ey a 


‘ . 2 
Pressure rise over a wavelength, Newtons/m 


ROLES 


Normalized pressure rise, Cine 


r1°.ITnstantaneous fluid flow, moeee 


Time average fluid flow, ne eee 
; - €Q> 

Normalized time average fluid flow, Po 
Time varying component of channel width, meters 
Equilibrium channel width, meters 

, ue 5) 
Fluid density, ke/m 
Time, sec 
Maxwell stress tensor, Neneene te 
Fluid velocity, m/sec 
Channel. width, meters 
System energy, watts 
System co-energy, watts 
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Notation Name and Unit of the Symbol 


X Transducer gap width, meters 
xX. tbe support thickness, meters 
y Spatial co-ordinate 

Z Spatial co-ordinate 

G) Frequency, 21ff, radians/sec 





CHAPTER I 
INTRODUCTION 

Peristalsis is the rhythmic motion of the walls of the 
alimentary canal, the ureter, and other hollow organs of the 
body, consisting of traveling wave contractions and dilations 
that move the contents of the tube onward, Outside the body, 
peristaltic pumping can be achieved by mechanically squeezing 
a fluid filled channel or tube in a traveling wave manner. 

The potential value of such a device must be apparent, 
In industry, the need often exists to pump fluids without 
having the fluid touch the pump apparatus, as in the case 
of blood and caustic chemicals, The peristaltic pump enables 
the insertion of a fluid carrying tube into the pump with 
no actual contact between the fluid and the ae mechanism, 
In addition, if such a device could be miniaturized, it might 
be inserted into the body to create artificial body functions 
when natural functions have fen 

All present peristaltic pumps rely on transverse punp 
tube deflections, which travel longitudinally, to create fluid 
Motion, \1'*?) Deflection is generally produced by a mech- 
anical roller which moves along the tube, squeezes it, and 
pushes the fluid onward, Electrical forces are Sneed only 


in that an electric moto} produces the force which causes 


the longitudinal mechanical motion, 





This investigation deals with a different class of electro- 
mechanical devices, In general, mechanical motion is purely 
transverse and electric forces are rectly related to pump 
tube eae ions.) Schematic representations of such devices 
are shown in Figure 1.1. In Figure 1.1.a, the channel walls 
form a continuum electromechanical transducer. The source 
of electrical force is a current sheet which is wound similar 
to the windings in an electric motor field and which produces 
a sinusoidal traveling wave. If the channel is an elastic tube, 
the continuum transducer can be approximated by arranging 
a large number of discrete magnetic segments on top of the 
tube, Their response to the traveling magnetic wave will 
cause tube deflections. 

Figure 1.1.6 represents a peristaltic pump which also 
moves purely transversely, but the sources of tube deflection 
are discrete transducers rather than an eect one chanical 
continuum. The schematic shows each transducer being one 
quarter wavelength, In order to produce a single wavelength 
of a traveling wave, as illustrated, each transducer excitation 
must be phase shifted ninety degrees from its predecessor, 

The number of discrete transducers per wavelength can be in- 
creased as desired, provided they are properly excited, In 
the limit of many transducers, the discrete pump approaches 


being a continuun, 





SCHEMATIC REPRESENTATION OF A TRAVELING WAVE PERISTALTIC PUMP 











x A 
Qe 
5 OOO@ 7 83 a a 2 
Current Sheet 
x } 
Transducer al. 
Actuator 
-* 
ee i, a pa ar) 
| | ibe 


| 
| 
sal = 
Cee oo a 


0 Transducer Exciters xX 






A>> 
A>>k 


Pressure = P(z) 


Velocity V(x) 


Corts) 


Pirure i.) 


tele. 





In both representations, Figures i.1.a and i.1.b, the 
electric force which causes tube deflection is a function 
of the gap width, X, between the actuator and the exciter. 
It will be shown that this relationship produces a mechanical 
instability which causes the discrete configuration of Figure 
1,1.b to operate in an on-off, that is, channel open-channel 
emt off, mode. 

The analytical model for the pump is based upon the physical 
situation most likely to be encountered where peristaltic 


(5) 


flow exists. Transverse tube dimensions are small compared 
to a wavelength of Chemereets nz wave, Transverse velocities 
and pressure gradients are small compared with longitudinal 
ones. The pressure in the tube is assumed to be instantaneously 
uniform over each cross section, The fluid density and 
viscosity are assumed to be constant. The Reynold's number 
is assumed to be low enovgh so that flow is viscous, Jn order 
to insure the validity of a viscous model, Reynold's numbers 
should be less than eo | 

The tube walls move in the transverse direction only. 
The traveling wave shape will be analyzed as both sinusoidal 


(7,8) 


- ands quare Pumping does not seem to be sensitive to the 
Shape of the wave, but specific values of flow rates are sensitive 
to the type of tube deflection, 


The pump developed for investigation in this work is 





shown in Figure 1.2. It consists of four transducers which 
are excited to produce a single traveling wavelength along the 
pump tube. The tube is inserted between’ the cates exciter 
and actuator. The actuator is supported by four springs, 

one at each corner, and by the fluid pressure in the tube 
passing underneath, The transducers have two exciter windings, 
One supplies direct bias current, The other supplies alter- 
nating current, The alternating current ee is a sine- 
cosine potentiometer, the output of which is amplified through 
amplidynes to produce the desired low Fer aie current values, 
The phase shift in excitation between each transducer is ninety 
degrees, 

The analysis of the pump is done step by step, starting 
with a single transducer and ending with the completed per- 
Peitic pmo, 

In chapter two, the inductance of a single transducer, 
the supporting spring constant, and the electric forces exerted 
on a transducer as it will appear in the pump are determined, 
Then the time variation sa single transducer and the set 
of four transducers in sequence is ascertained for various 
direct and alternating current excitations. The frequency 
of alternating current is arbitrary except that the Reynold's 
number is directly proportional to frequency. Lower Reynold's 
numbers result from lower frequencies, The lowest possible 


Reynold's number is desired to have viscous flow forces dominate 


13. 
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PERISTALTIC PUMP PARAMETERS 
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inertial flow forces, 

Once the transducers are analyzed individually, the per- 
istaltic pump of which they are a part is analyzed as a whole 
in chapter three, First, an empirical expression for the 
pressure force exerted by the pump tube is obtained, Then 
a theoretical expression describing transducer actuator motion 
is determined, Finally, based on known actuator deflections, 
the pressure-flow characteristics of the pump are determined 
for a square traveling wave and a sinusoidal traveling wave 
description of tube deflections, The experimental operating 
points for the pump are compared to these characteristics, 

Throughout chapters two and three, analytic work is 
immediately followed by test results, This is done because 
experimental results from some sections are used in following 


ieyolc Work. 
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CHAPTER IL 
MAGNETIC TRANSDUCER ANALYSIS 
2.1 Magnetic Transducer Inductance 
Lote Theoretical Inductance 


The theorétical inductance of a single transducer is det- 


(9) 


Apply Ampere's law to the transducer as pictured in Figure 2.1la. 


ermined using lumped parameter electromechanical techniques, 


The permeability of the exciter and actuator are assumed to be 
infinite, This requires that the magnetic field intensity, H, 
must be zero since the flux density, B, is finite and B =A HL, 
The only magnetic field is in the gaps between the exciter and 


the actuator, 
\,V-X H+ nda = \, J + nda 21a 
Apply Stokes’ theorem to obtain 
§J.H + dl = §.d + nda : Bails 


The contours for these integrals are indicated as dotted 


lines on Figure 2.1a. Integration yields 


H, Ati xA=N 1 Ze 


2 


Taking a surface that encloses the transducer actuator in 


gee 
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Peeure 2,16, Knowing that VY« B = 0, integrating VY: B over the 


enclosed volume, and applying the divergence theorem yields 
§.B + nda = 0 2S 


In terms of the magnetic field penetrating the enclosing surface, 


as becomes 
Zu, Hy S| ie - de He, b a> =20 24 


Solve equations 2.2 and 2.4 simultaneously, 





b NI 
oy” 3a | heee 203 
La N I 
to “Zat+b Xx 2.46 
The total flux through the center leg is 
@ =4, H, bd 25 
_ 2Ko iN] Gael ey ve > 8 


X (2a + b) 


In the absence of leakage, the flux links the N turn winding 


N times. The flux linkage is 


Os 





N ® 2S 


ed 
1 


2M wrabd 


Ue ee 2.10 
X (2a + b) 
= iGO Pol 
21 RS 
L(x) = “0 2.12 


eres 


X (2a + b) 


2.1.2 Experimentally Measured Inductance 


The transducer exciter has a primary winding of 500 turns 
and a secondary winding of 26 turns on the center leg. A com- 
parison is to be made between the theoretical transducer inductance 
and the measured primary and secondary inductances, Because the 
secondary Snductance winding has only 26 turns, it must be nor- 
Malized in proportion to the turns ratio squared if it is to be 
made comparable to the primary inductance on the same plot of 


inductance versus gap width, 


_, 500 \2 
L(X) normalized — ( 26 ) LX) . measured 2.15a 
= 370 L(X) __. | 2.13b 


The theoretical transducer inductance is evaluated from 
equation 2.12 by substituting the transducer parameters from 
Reese 1.3. 
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?)y; ‘-L} 
r nex ime as x 10) . 
LX), ra 4 ee 14 


Plots of transducer inductances versus gap width are shown 
in Figures 2.3 and 2.4, Further comparative insight can be 
gained by taking the natural. logarithm of equation 2.14 and 


differentiating with respect to in X, 


in L(X), ea lay aot oe 107" - ln xX | 2 ella: 
din LX), 2 24 2,15b 
Clavie ak 


A logarithmic plot of theoretical transducer inductance 
versus gap width is a straight line with slope of minus one, 
Figure 2.5 shows eee ena plots of theoretical mepraimary, 
and normalized secondary inductances, 

The measured values of primary and secondary inductance vary 
Peeiderably from the theoretical values, Tyeeeetion of Figure 
2.11 shows that electric forces calculated using the theoretical 
transducer inductance are closer to the actual value of measured 
electric force than electric forces calculated using measured 
inductance values, The range of gap width for these calculations 
ais from ,003 to 00! meters, In the limiting cases of X = 0 and 
X + infinity, the measured inductance values are more accurate 
than theoretical ones, Theoretically, as X approaches infinity, 


the inductance approaches zero, And, as X approaches zero, the 


Zhe 





2°? GANSTA 


SHONVLORGNI UAOACSNVEL TVOTIZYORHL UNV GEYASVEN ATIVINGWIUSEXe 


61000° , oD 

22000° ¢1600° 

€Z000° : 08500° 

6Z000° ¢6Z00° < 


9£000° €2T00° 
$S000° 977000° 


Z6000° 


4. | us us ; 
; d 
(X)T Cn | (X)T (x) y 


WONVIONGNI | sONVIONGNT FONVIONGNI 
SJONVIOOANI XEVGNOOGS CTUNSVAN CUIASVAW 
TWOLLHUOREL CAZTTWRUION RUVANOOAS xavd HIGIM dv) 








Z00° 


(stoqzow) x ‘uyptm dey 


900° GO0° 700° €00° Z00° Loo © 





*ex)4 


HLGIM dVO SQSYAA SHONVIONGNI.ugoOnasnvas TWOLLHUOUHL -GNV “RUVGNOOES ‘AUVWTYd HO Lotd 


€°0 


cqusy). (X)JT feourqzoupuy 


(so 


Pace 3 





EXPANDED PLOT OF PRIMARY AND THEORETICAL TRANSDUCER INDUCTANCES 
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Inductance, L(X) (henries) 
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inductance becomes infinite, Practicelly, neither of theso the- 
@revical limits is acevrate, In the real infinite limit, the 
gap width is so great that the magnetic field lines find a path 
of less reluctance by going directly from the center exciter leg 
to the outside legs through the air. Even in the case of less 
than infinite gap width, some field lines Grave lida rec uly om 
leg to leg altering the transducer inductance. In the zero limit, 
when no air gap exists between the actuator and the exciter, the 
actual inductance is not infinite because the permeability of | 
mae actuator and exciter metal is not infinite, as was assumed 
in the theoretical model, but finite. 

Fringing fields and non-infinite transducer permeabilities 
affect the measured inductance values for all gap widths. Field 
' losses can account for the difference between primary and sec- 
ondary inductances, The primary exciter winding is wound directly 
to the exciter center leg. The secondary winding is wound over the 
primary. The tendency would be for greater losses to be experienced 
in the secondary, resulting in a eeeee inductance. Pigures2 .3 
shows this to be true, 

Additional errors could have come from the impedamce bridge 
used for inductance measurements, 
2.2 Actuator Support Spring Constant 

Each transducer actuator is supported by four springs, The 
value of the spring constant, K, for these springs was desired, 


Weights were placed on top of each actuator, The resulting 


2a, 
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PLOT OF SPRING FORCE ACTING ON THE TRANSDUCER ACTUATOR VERSUS SPRING 
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change in gap width due to these weights, which is the spring 
deflection, was measured with micrometers. 

The values of weights and spring Perlecmicn are plotted 
On Figure 2.7. The spring constant for the four springs acting 
in parallel is the slope of the curve, Thus, the spring force, 


f_, is expressed as k (Xp -~ X) where X, is the undeflected gap 


@) 
width. Figure 2.7 shows the spring constant to be 2860 newtons / 


meter, 


meovuillectric Forces 


el Electric Force Determination Using Energy Methods 


(10) 





The magnetic transducer is represented as a two terminal 
circuit in Figure 2,8a. The total energy stored in the system 


is Wa The power or rate of change of energy is 


oe ee f. = 2.16 
at 
eocne dx 
~ 7 at te dt 217 
dW =iddA -f. ax Danes 
m = : 


it is desired to have current as the variable, sco 2 co- 


energy is defined as 


dWt=-dW +ind ys 
m m 
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b, Paths Of Integration For Evaluation Of Magnetic Transducer 
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pillgstitute equation 2.19 in 2.18. 


dwt = \di+f_ dx Ze 
mA € 
Integrate the co-energy along the path shown in Figure 2,5b, 
The integral is independent of the path taken, Along path G) : 
Gurrent is zero so the electric force is zero, Along path (2) : 


tae displacement is at its final value, X. 


x ( 
aoe : 
in §, Adi 2,21 
C=0 y=\¥ 


Using the value of A from equation 2,10 and the definition 


emeex) from 2,12, 


Zz 
(age Sn nee 
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251) D238 
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In differential form, dw can be represented as 


t ye 
amir Gav: a ax RES 
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i Git. ele 


Compare this equation with equation 2.20, 





d t 
f, = es 2.26 
AG 


rom equation 2.22, the electric force is 


2 
¢ --AWaba y See 
1 (2a + b) 
: qd) ie 2.28 


2e30e2 Electric Force Determination Using Maxwell 
cay 


Stress Tensor 

The electric force on the transducer necreitor in thee Greer on 
is desired, This is evaluated using the definition of the stress 
tensor in Figure 2.9. Since the only magnetic fields are in the 


x direction, only an x component of the stress tensor exists, 


ne ME 
T= 3A, CHT + HS + HY ) DPS 
fp}. Div nda 2.30 
ee 2 2 
Se eC 2 se) Gl a URIS 1) Gl) 265 
Substitute equations 2.5 and 2.6 for H, and Hy. 
Mk Waba 
f = 0 j 2 ZEe 
ee ere 
X~ (2a + b) 


Piohem 16 identical to equation 2.2/7. 
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MAXWELL STRESS TENSOR REPRESENTATIO! OF THE TRANSDUCER ACTUATOR 


ELECTRIC FORCE 
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2e3e3 Experimental Evaluation Of Electric Force 
For any applied current to the exciter windings, the actuator 
will assume an equilibrium position such that the springs are 
compressed and exert an upward spring force which is balanced 
by the applied electric force and the weight of the actuator 
acting Besnard 612 
The actuator weight is balanced by an equilibrium spring 


force, 
Me a ( Ko a xX, ) Ee 


The electric force will then be balanced by the additional 


compression of the springs. 


ye = ( x, =e K ) 2.34 
This force is the experimentally measured electric force. 
The electric force can be determined in two other ways. 


iimcoretically it is 


t= TOKO meee 
dx ie 


from equation 2,28 where a LOO is the slope of Figure 2.4. 


~ 


By use of measured inductances and the slope of the measured 


inductance curve, measured electric forces are 


34 ¢ 
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Over the range of gap widths measured, the theoretical electric 
force most accurately predicts Uhevectual electric orce, t mt 
The electric force calculated from the measured inductance is 
considerably in error. Jt was evaluated using the slope of the 
measured inductance curve, Figure 2.44, from section 2.1.2. An 
error might have been expected for the ee previously 
mentioned in section 2.1.2. 

2.4 Magnetic Transducer Force Eavili.briun 
Newton's law applied to the transducer actuator is 
M ax 2 8S = 1 ot ~ Mg Desi 
ate Ai spring electric 

The position that the transducer actuator takes in steady 

state equilibrium represents a balance between the upward spring 


(13) 


force and the downward electric and weight forces. 


F(X) 


i 


f + ie = Mg DNs 
S e 


2 
K(x, -x) + Ho WDE = ies D9 


ie (2a + b) 


iH 


N 


2 2 
N~abd I 
K (XX) +k (h,- x) 4 HM, Wabd Tove 2,40 
x* (2a+b) 


The weight of the actuator cancels K (Xp -X)). 
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F(X) = K (x, - x) + Ho Ne aba 1 2.44 
Xo (2a) 

Pagtre 2.15 1s a plot ef these tae functions eo: A.) ibe 
points where they intersect are the equilibrium operating points 
or the equilibriwn positions assumed by the actuator. A separate 
plot will be obtained for any value of current. The points, Xx 
are unstable equilibrium points, A change in X from this valve 
causes a force imbalance which tends to cause X to change even 
matcher, The points, Kos are stable, Any change in X from Xe 
causes a force imbalance which tends to return the actuator back 
18) X oe The actuator always takes a stable position, 

As current is increased beyond 0.677 amperes, the values 
of electric force are always greater than spring force, There 
eee. NO equilibriun positions for Chevsoula vor lectures croc 
always overcomes the springs and closes off the gap, The point 


of tangency, X,,, is the equilibrium position for the actuator 


di 
when current is equal to 0.677 aimenoe: Any perturbation acting 
to increase X causes the actuator to return to Kove Behe hens 
perturbation which acts to decrease X will cause the actuator 
memclose off to X = 0. 

The accuracy of Figure 2.13 can be checked, Each stable 
equilibrium point, X os on Figure 2.13 should sennccuora to the 


measured gap width from figure 2.10 for the same value of current, 


Pecamination shows this to be true, 
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PLOT OF SPRING AND ELECTRIC FORCES VERSUS GAP WIDTH FOR VARIOUS CURRENTS 
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Time Variation Of Transducer Actvator Position For Sinusoidal 
Poplaed Current | 

Plots of actuator position versus time are obtained using 
the current-gap width correspondence from Figure 2,10. Figures 
Z.i4a, 2.i5a, 2.16a, and 2,17a are plots of time varying applied 
current and time varying head position for various combinations 
of direct bias current and alternating current, The gap width 
"closes off" at X = .00173 meters because the tube support spacer, 
as shown in Figure 1.2c, prevents it from going to zero, When 
the pump tube is inserted, it will be squeezed completely shut 
inecnis condition, 

It was desired to operate the pump less than fully occluded, 
In concept this could have been accomplished by limiting the exciter 
currents such that the actuator never fully closed off the channel. 
Such operation is illustrated in Figure 2.17a. Under these conditions, 
the channel has a steady state width and a time variation in 
width around the steady state, Schematically this is shown fe 
Figure 1.1 when R>r. Practically, operation of the pump as 
described is very difficult. During operation, variation in 
the characteristics of the individual transducers from each 
other cause one to operate properly while another closes off, 
Figure 2.17c shows actuator three operating almost paneccey 
maile the other three are closing off, The transducers are so 
mistable that only slight porturbations ere one to make 


the gap close off completely. 


Ly 





Zo cvercome this difivculty, a channel dimiter was pilaced 
astride the pump tube, Actuator motion is restricted because the 
actuator contacts the rigid limiter as it tries to close off, 

The pump transducers can be driven as hard as desired, but peristalsis 
takes place only on that section of the tube above the limiter, 

The thickness of the limiter can be varied to obtain the desired 
values of Rand r, Figure 2,16a illustrates the time variation 

of actuator motion with the limiter in place and Figure 2.16c 

shows the Sanborn recordings of actuator motion under these 
circumstances, 

Figures 2.14b, 2.15b, 2.16b, and 2.17b are representations 
of the positions of the four actuators for various times covering 
one complete cycle, Each actuator is phased ninety degrees from 
the preceeding one, The resvlt is a traveling wave to the right 
which can be approximated as a square wave or as sinusoidal, 

Gamalily, Figures 2.i4e; 2.1 5e,02e1 cewmenan 247 oem tne 
which have already been discussed, are Sanborn recorder plots 
OF actuator position versus time for all four transducer actuators, 
The plots are from the output of photo electric cells which were 
arranged to sense the motions of the actuators, Determination 
of the exact magnitudes of gap widths is not possible from the 
recordings because the recording amplitudes are relative and 
a function of the sensitivity setting, However, the recordings 
do illustrate the Paes hice between the actuators which creates 


the traveling wave, The phase shift is not ninety degrees between 


WO. 





each actuator, The outputs of the power supplies for the exciter 
currents were not quite properly phased, The resulting traveling 
wave was mete icilesa close pause Biggs cies pause close 
BICSE. occ eceveoo 

A comparison between the predicted actuator time varying 
position (e.g. Figure 2.15a) and one cycle of an actual measured 
time variation (e.g. Figure 2.15c) shows that the actual motion 
conforms quite well to predicted motion, The operation of the 
transducer is effectively open-closed, The transducer is so 
unstable that it acts almost like an on-cff relay, This leads to 
the mension that an elaborate multi-phase sinusoidal power 
supply is not necessary for efficient operation of the pump, 
A simple on-off direct current power supply for each transducer 
which is actuated ninety degrees ahead of its predecessor is 
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CHAPTER fA 
PERISTALTIC PUMP ANALYSTS 

Analysis of the electrical and mechanical properties of 
individual transducers has been completed in tne previous chapter, ° 
This chapter is concerned with the analysis of the peristaltic 
pump as a whole. The fluid filled tube running between each 
transducer actuator and exciter exerts an upward pressure force 
on the actuator that must be determined and entered into the 
force balance on the actuator. The magnitude of this upward 
pressure force relative to the electric force, which acts to 
pull the actuator downward and close off the flow channel, must 
also be determined for the following reasons. If the pressure 
force is comparable to or greater than the electric force, the 
periodic motions of the transducer actuators (i.e, tube deflections) 
are not straight forward to analyze. Fluid mechanical effects 
between transducers must be accounted for, If the electric force 
dominates, the tube deflections have the same een as the 
applied magnetic traveling wave and approximate analysis of tube 
deflections is simple. 
jeeecump Fluid Pressure Forces 

3.1.1 Theoretical Description Of The Pressure Force 
Exerted On The Transducer Actuator By The Pump 
Suction Head 
Pie stearic 2 1uid Seceeene exerted on the tube wall and the 


transducer actuator is 2a function of the suction head, H Tnertial 


® 
S 
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or fluid velocity effects on the fluid pressure will be neglected 


in the analysis, 


eet oes 3.1 
The total fluid pressure force acting on each actuator is 


i alt 
p Ww Cae 


= Pewld, 3,3 


where w and | are dimensions from Figure 1.2, 


wolve equation 3.3 for wl, 
So 


wl is the variable area over which the fluid pressure acts 
Seerne actuator, It is a function of the gan width and will 
be derived experimentally in the next section, The emperical 
expression for wl would be used in déscribing the motion of the 
peristaltic pump, section 3.2, under operating conditions where 
the fluid pressure is comparable to the electric force, 

3.1.2 Experimental Evaluation Of The Fluid Pressure 
Force 
Figure 3.2 is a plot of the total force due to springs and 
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Figure 3.2 





fluid pressure acting on each actuator versus channel width, 

Since the spring Moraes is known, it is subtracted from the total 
force, leaving the force due to the fluid mes ceuie alone, The 
fiuid pressure curves yield an empirical expression for the fluid , 


pressure force as a function of suction head and channel width, 


a Se eG He - 750 (X - xX.) Bao 
where A. and X are in meters and f is newtons, 

Substitution of this empirical expression in equation 3,4 
yields 


1 2a He SO = xX.) 


Vl 2 Se ee Sie 
Pek, 


Comvarison Of Relative Magnitudes Of Forces Acting 





On The Transducer Actuator. 

The forces acting on the actuator vary with gap width, suction 
head, and exciter current. The fluid and spring forces act upward, 
while the weight of the actuator and the electric force act 
downward, During operation the gap width varies from about ,004 
to .0017 meters. Table 3.1 shows values of forces for different 
gap widths, suction heads, and exciter currents. 

The figures in Table 3.1 do not clearly indicate what happens 
when the pump operates. As current to the exciter reaches its 


peak value, the actvator begins to close due to a net imbalance 
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abpronate Force Values For voeus Gap Widths 
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of forces acting downward, As the actuator closes, the electric 
foree increases at a rate greater than the spring and pressure 

forces combined, The actuator is unstable. Once it starts to 

close, the imbalance of forces increases and it will close completely. 
Further, as the channel pinches closed, all fluid is forced out 

and the fluid pressure drops to zero, Observation of operation 

of the pump shows that the fluid pressure and fluid dynamics 

do effect the pump, but that the arose cone does dominate 

and the periodic motion of the actuators has the same frequency 

me the applied exciter current, 

Dee theoretical BeciereLaeelein Of Deietal veer anebane enon 115) 


The balance of forces acting on the actuatcrs in Figure 1.2 


yields onl 





P (wl) + K (Xo - X) = Me 4 - Bay 


In static equilibrium, the weight is balanced by an equilibrium 


Spring force, 
Thus 


P (wil) + Kk (X, = X) = we oe 


Rearrange equation 3.9 and divide by wl. 


gees OE 


a 
wl wl 


SIRS: 


Substitute for f, from equation 2.27 and for wl from equation Sree 


NO 


pels 2 Oe 
7 ee om Tk Gt, - 0] ae 


>< 


Consider the tube wall to be a continuously deforming surface 
excited by a traveling wave. That is, the discrete motions of the 
transducer actuators are approximated as an appropriate continuous 
traveling wave, 

The tube wall deflection travels to the right in Paeure 3.5 
with a constant velocity, c, which is a function of the frequency 


and wavelencth of the traveling wave, 
G2. 





CROSS SECTIONAL VIEW OF THE PERISTALTIC PUMP IN THE FIXED AND 


MOVING FRAMES 





eee 
i 02 at Ce 
P' = P'(z) 
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V Vix) 
r = r(z') 
Figure 3.3 
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c=)\f 3 Baie 


where ) = 41 from Figure 1.2 
The tube wall itself is fixed in longitudinal position and 
moves only transversely. The tube wall position in the wave or 


moving frame of reference is 
Oe REC al | bos 


In the wave frame, the traveling wave is a periodic stationary 
deflection, By assuming steady state flow conditions and ignoring 


inertia, the Navier-Stokes fluid flow equation is 


VP = Ta Aan 3.14 


The model for the tube or channel is quasi one dimensional, 
The pressure is independent of transverse position. The velocity 
is purely in the longitudinal direction, but varies with transverse 


position, 








| 
dP’ .o 3.15 
‘ebay 
d pl _ cans Bae 
dz' Ah 5 ‘ 
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In the wave frame of reference, the volume rate of flow 
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per unit depth in the y direction is 


X 
Q'= ii ch Dollie = 
5 (2 


The fluid velocity is determined by integrating equation 


Bato and applying the boundary conditions 


I 


ve (x =X) =v_ (x =X) -c = <c | Srcike: 


I 
© 
w 
I 
t 

Q 


Vv, (x 3.19 


ime result is 


_ 


er oa Ge GS) Se 3.420 


Substitute equation 3.20 in equation 3.17 and integrate 


to determine the expression for flow rate, 


1 @P 3 


Pe ; | 
in correo) Boel 


Rearrange and solve for the pressure gradient, 


ar = - TAs (0 + eX) Bees 
Az x? , 


Now solve for the time average flow rate in the fixed frame 
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with the ultimate object of obtaining an expression for the flow 
rate in the moving frame for substitution in equation 3.22, 


| x 
<a> =<] v, ox > 3.23 


x 
-<f (vi +c) dx » 
0 


ue farst part of the integra lgac.O.. 


phase velocity times R, the time average value cf X, 


Bee 


The second is the 


3.25 


<G> 10 ok 


Rearrange. 


NaAe 


OY = (G7 — cx 


which is independent of longitudinal position, Substitute Q' 
= Or, . 


in equation 3.22 and note that c (X - R) 
BAY 


dP 12 
ee a i6Ge tee) 


The expression for pressure from equation 3.10 can be dif- 


: d d 
ferentiated with respect to z, since es ee to yield 
K (X, - X 
1244 (<Q> + er) = S| fe se ital 3,28 
. dz wl wk 


Xx 
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For a given periodic driving current, equation 3.28 represents 
a non-linear equation which can be solved for X. Only one boundary 
condition is required, The response desired is periodic with the 
same wavelength as the driving current. The required solution 
for X is found by adjusting the boundary condition at one position 
such that the deflection returns to the same value after one 
mavelength, that is , X(z‘) = X(z" + 4), 


Determination Of Pressure-Flow Characteristics For A Known 





imeansducer Actuator Deflection 
3.3.1 Description Of Peristaltic Traveling Wave Deflections 

The peristaltic pump developed for this thesis, as stated 
in section 3.1.3, operates so that electric forces dominate, The 
traveling wave frequency is the same as the anplied electric exciter 
frequency. The transducers come very close to operating as open- 
close relays. Thus, an approximate description of the tube deflection 
is known. The non-linear equation 3.28 need not be solved to 
determine X, Rather, the SNe ss uneee ler characteristic for the 
pump can be determined by integrating equation3.27 under the assump- 
tion that all quantities inside the integral are known functions 
ot Ze 

Figure 3.4 is a Sanborn Recording of the tube deflections 
during pumping. This figure can be compared to moeure Zell ger. fal 
similar trace except that no fluid filled tube is in the transducers, 
The erratic traces on Figure 3.4 are due to sensitivity adjustment 


difficulties on the recorder, not fluid action, The presence of 
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fluid does not significantly change the traveling wave. 


, d ap | 
AP = [ pare dz Chere 
d a dz 
d <Q>+ cr 
- 2a {| <a dz Delo 
O xX . 


The peristaltic traveling wave can be approximated as a 
square wave or aS a sine wave, Both approximations are made 


in succeeding sections, 


3.3.2 Square Traveling Wave Approximation To 
; ele) 
Perastaltic Pumpin 

fhe peristaltic traveling wave can be approximated as a 
square wave. This is illustrated in the laboratory and wave 
frames in Figure 3.5. 


The pressure rise over a wavelength for a square wave is 


determined geometrically using equation 3.30. 
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OQUARE WAVE PERISTALTIC PUMPING 





a. Square wave deflection in the moving frame, The wall velocity 


is -c and the volume flow rate is independent of z'. 


0 


b. Square wave deflection in the fixed frame, The increasing 





volume in the half wavelength between 2 and } necessitates 
a net flow into that volume, The velocity profiles would yield 


a net flow out of the decreasing volume between 0 and Agee 


Figure BG 
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Substitute these values in equation 3.31. 


Eps = CQ? +h :e Eager | 
_ (1 +h)? (1 +h)? 

The normalized pressure rise as a function of normalized 
time average flow rate is plotted in Figures 3.6, 3.7, 3.8, and 
3.9.for various values of h. As h approaches one, the channel 
as cutoff and the pressure rise becomes ee Such a positive 
displacement pump, operating so Sis channel is cut off while 
pumping, can theoretically pump against an infinite head, 

When h is zero, there is no channel wall deflection. If 
a flow rate existed in the channel, this Ses would give the 
pressure rise up the channel for a given flow rate with no wall 
deflection, However, in the specific case of the peristaltic 
pump where all flow is due to wall motion, the flow rate and 
pressure rise would also be zero, 

The intercept points for the normalized pressure-flow curves 
are determined by evaluating the normalized pressure rise for 
zero flow rate and the normalized time average flow rate for 
zero pressure rise, Intercept values for operating suction heads 


and normalized channel width variations are shown in Table 3.2. 
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oquare Wave Pressure-rlow Intercepts 





3.3.3 Sinusoidal Traveling Wave Approximation To 


Peristaltic Pum meh 


The peristaltic traveling wave can also be approximated as 
sinusoidal, This is illustrated in Figure 3.6. The pressure 
rise over a wavelength for a sinusoidal traveling wave is 


determined by integration of equation 3,30. 


AP _ 124 A <Q> + cr cos(vt -kz) ( “k aie) B28 
A “ o [R + r cos (wt - ka)} 7: 


When Sey this integral becomes 


= CPST 3ckr~ - £Q > (2R° + | 3.639 
ate 


Normalize equation 3.39 using the same normalization as in 


the previous section, Then 
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SINUSOIDAL WAVE PERISTALTIC PUMPING 





a. Sinusoidal wave deflection in the moving frame. The wall 


velocity is -c and the volume flow rate is independent of 2', 





0 d ees 


b, Sinusoidal wave deflection in the fixed frame, The increasing 


volume in the half wavelength between )/2 and i necessitates 
a net flow into that volume, The velocity profiles would yield 


a net flow out of the decreasing volume between 0 and 3/2, 


Figure 3.6 
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The discussion in the previous sevtion is applicable to 
the sinusoidal approximation of the peristaltic traveling wave. 
As from the previous section, the intercept points for Figures 


3.7, 3.8, 3.9, and 3.10 are below. Specific values of the points 


are in Table 3.3, 


ge 2D _ 3.41 
2 | 
2a 





Table 3.3 


Sinusoidal Wave Pressure-Flow Intercepts 
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3.3.4 Experimental Pump Operation 


The pump was operated peristaltically with the same suction 


head and normalized channel width combinations for which pressure- 
Hlow characteristics were determined in sections 3.3.2 and 3.3.3. 
Values of nurs flow rates and pressure heads were measured and 

are recorded in Appendix I, The operating points are plotted 


on the appropriate pressure-flow curve, Figure 3.7; 3.0, 3.9; or 3.00. 
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NORMALIZED PRESSURE-FLOW CHARACTERISTIC FOR PERISTALTIC PUMPING 
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In Appendix I, the measured flow rate was divided by the 
channel width prior to normalization to obtain flow per unit 
channel width, The pump pressure differential was determined 
by measuring a static height difference between the pump suction 
and discharge headers, The normalized flow and pressure differ- 
ential are determined using the same definitions as in sections 
343.2 and 3.3.3 and applying the measured parameters of frequency, 
wavelength, channel width, and phase velocity. 

3.3.9 Cut Off Or Occluded Pumping — 

As was pointed out in section 3.3.2, when h =i, the channel 
is cut off during pumping and theoretically the pump can purap 
against an infinite head, In this case, a pressure-flow char- 
acteristic is meaningless, 

Theoretically, the,entire contents of the pump are pushed 
forward at the phase velocity of the applied traveling wave during 


occluded pumping. For a unit tube width, 


Using the same definition for normalized time average flow 


moe 2S in Section 3.5.2; 





<Q>. 
<Gy = Re 3645 
es GaP ere 3.46 
<e? theoretical 
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<Q> can range *»om zero, no flo;, to one, perfect occlusion, 
The experimental normalized time average flow rate in Appendix II 
ranged from 0.26 to 0.45, Considering the sources of error, 
which will be mentioned in section 3.4, the most important to 
inefficient occluded pumping are improper time phasing between 
actuators and the fact that the pump consists of only one wavelength, 
Both of these items bring about physical conditions under which 
the entire contents of the tube under each actuator are not pushed 
forward as the actuator squeezes downward, 

If the pump were operated so that each transducer was on or 
oir, Almost ideal occlusion could be ebtainea aa each actuator 
was held down until the succeeding one was closed, Under these 
circumstances, the fluid has no choice except to move onward. 

For the specific pump ir this thesis and under the type 
of operation finally used where maximum current was applied for 
all pumping, peristaltic and occluded, occlusion is the most 
efficient form of pumping. For a given power iipuL, Occlus iON 
gives the greatest flow rate. Unless less than fully occluded 
pumping is desired, the most effective pump would be built to 
operate fully occluded, 


a ae Se em ee oe ce 


The appropriate Reynold's number describing flow in the 


(19) 


peristaltic pump has been shown to be 
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The Reynold's number eee between |,31 and 8.65 for the 
operating points of the pump, Exact ae are shown in Appendices 
I and II, In each case, inertial forces nre greater than viscous 
forces, In spite of this, the viscous malel developed yields 
reasonably accurate results. | 
3.4 Sources of Error 

Hits Ge Ve Ulaeero attribute the error between theoretical. 
pressure-flow characteristics and experimental operating points 
to any single cause, The inertial nature of the pump, which is 
not accounted for in the theoretical mode), could cause some 
error, The mechanical aspects of pumpin;, could have ‘produced 
errors, The rubber pump tube was somewhat elastic which could 
have altered flow, The tube was eliptical in shape rather than 
a square Saararel as the theory assumes. The model assumes the 
vertical walls of the tube are always verlical and perfectly 
elastic to up and down motions, while the channel. is of constant 
width. The real pump has elastic vertical eee eels bow out, 
Channel width varies as the pump actuator moves up and down, As 
was pointed out in section 2.5, perfect ninety degree phase shifting 
between transducers was not obtained. The uneven actuator motion 
could have affected flow, The pump pressure head, from which 
AP results, is on the order of one millimeter, Accurate mea- 
surement of these differentials was difficult. The pump has 
macher lone sections of aeons tube on Lhe inlet and outlet, 


The flow resistance and surface tension «tfects at the ends of 


Si, 





these tubes could alter pressure-flow relations, 

Considering the sources of error, it is surprising that the 
operating points are as close to the theoretical pressure-flow 
curves as they are, But, the real value of the project is not 
the relative agreement between theory and experiment, Its value 
is in the physical realization that longitudinal peristaltic 
type flow can be produced by purely vertical transducer motion, 


No longitudinal mechanical motion is involved, 


82 





CHAPTER TV 
CONCLUSIONS 

The experimental work done in this thesis has shown that 
the observed peristaltic pump operating points lie reasonably 
close to the theoretical pressure-flow characteristics of the 
pump, It has not been shown that the eee eros alan Gcharacteristies 
accurately predict pump operation, Hirst, very few overating 
points were measured, No real trend or correlation between 
measured and theoretical operation could be observed, Second, 
the large number of error sources make it difficult to differentiate 
flow variations due to theoretical prediction from flow variations 
due to experimental error. The exception e this is occluded 
pumping. Complete occlusion theory predicts a normalized flow 
rate of one. Flow rates less than one are due to inefficient 
or improper experimental operation of the pump and are readily 
recognizable, 

The thesis and experimental work have been of anne. The 
major reason for the effort was to show that peristaltic pumping 
could be produced by a principle which had never been tried. 
Longitudinal flow was produced where no longitudinal mechanical 
motion was involved, The importance of this fact cannot be over 
emphasized. The longitudinally traveling veristaltic wave was 
produced by properly phasing the transverse motion of discrete 
electromechanical transducers, Jf a peristaltic pump is ever 
to be miniaturized and used as a replacement for an organ in the 
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body, the use of some sort of mechanical roller is not practical, 
But, such a pump could be designed, using the principles of this 
thesis, which would be easy to miniaturizo, mechanically simple, 
flexible, light weight, compact, and reliable. For instance, 

if very small displacements are desired, the transverse transducer 
motion could be produced by fer Cais ene ances of ceramic 
Seystals, The excitation for a piezoelectric crystal is a Voltage 
applied across it. 

Rete: new principle developed in this work is the direct 
relationship between the applied electric force and the transverse 
motion of the transducer actuator, In the case of the peristaltic 
pump developed in chapters two and three, the magnitude of the 
. electric force is inversely proportional to the gap width squared, 
In future So vane devices, it might be desirable to utilize 
the electric force or the magnetic field directly to produce 
fluid motion, This has possibilities in pumping magnetic liquids 
directly or in pumping by means of magnetic liquids, A traveling 
magnetic wave might be made to pull magnetic liquid along just 
as the traveling field of a motor pulls the armature along. Both 
of these principles are of value for future experimental work, 

The specific pump developed could be improved in the following 
ways. For most efficient operation, it should pump by totally 
cutting off the flow channel in its transverse motions, This 
causes the greatest volume flow rate for power input, However, 


total ceclusion covld be disadvantareous. Blood is very sensitive 
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to pumping, The effect of being pinched in the flow channel could 
damage blood so that it might not be usable, The pump should 

be improved for pumping both fully occluded and peristaltically. 
Phasing and phase sequencing need improvement. Proper application 
of excitation could produce rhythnie heart-beat-like pumping 

mauner than continuous smooth flow, sit that were desired, Since 
the transducers act almost like on-off relays, a complex sinusoidal 
excitation is not needed, More efficient, simpler excitation, 
which ae just as effective, could be provided by a DC on-off 
power Eee 

One of the most serious defects in the pump is the low volume 
flow rate, The pump tube used was a long, thin balloon, It was 
difficult to find an elastic tube which would rit in the gap 
between the actuator and the excitemiand which did not exhibit 
elastic forces opposing gap close off greater than the electric 
forces applied, Ideally, a long, very flexible, broad tube is needed 
which would fill more of the channel width and have negligible 
vertical wall rigidity. Total channel "elasticity" should be due 
to the suction height of the fluid in the channel, 

The type of electromechanical peristaltic pump investigated 
in this thesis is rugged, inexpensive, and simple, It was built 
using common transformers, springs, plastic, and balloons, The 
pump operated consistently and dependably, These features , and 


those already mentioned, make it worth additional developient, 
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